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Abstract 

The Navier-Stokes equation in the bidimensional torus is consid- 
ered, with initial velocity in the Besov spaces B P r +2 r and forcing 
term in L r (0, T; B~*) for suitable indices s, r,p, q. Results of local ex- 
istence and uniqueness are proven in the case — 1 < — s + 2 — | < 
and of global existence in the case — s < — s + 2— ^ < 0. 
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1 Introduction 

Analysis of classical or weak solutions to the Navier-Stokes equation in a two- 
dimensional bounded domain has been widely considered in the literature. 
A lot of work has been done on solutions with finite energy, because of 
their physical meaning. However, infinite energy solutions are important 
for other reasons. On the one hand, analysis of self-similar solutions in 
the whole space involves velocity fields with infinite energy because of the 
no rapid enough decreasing to infinity (see, e.g., [Cat IGPj and references 
therein). On the other hand, invariant (equilibrium) Gibbs measures are 
known for two-dimensional hydrodynamics in bounded domains (see, e.g., 



[ARH-K1 IAH-K] for the deterministic Euler equation and [DPD1 lAFbj for a 
stochastic Navier-Stokes equation, whereas |ACj deals with both the viscous 
and inviscous problems). The velocity fields with finite energy are negligible 
with respect to these invariant measures. This makes interesting the analysis 
of infinite energy solutions in a bidimensional bounded spatial domain. 

The aim of this paper is to investigate the Navier-Stokes evolution prob- 
lem in a bounded domain D of M 2 with initial velocity and forcing term of 
low regularity. The force will be integrable in time to some power with 
value in some Besov space B^ q of negative order a (for the space variable) 
and the initial velocity will belong to some other Besov space of negative 
order. We point out the importance of the force; for instance, this allows to 
consider problems related to the Navier-Stokes equation with a stochastic 
forcing term. (See, e.g., [DPDj for the study of auxiliary (deterministic) 
equations of Navier-Stokes type arising from a stochastic Navier-Stokes 
equation. There the forcing term has space regularity of a negative order 
Besov space. We also refer to the bibliography of [DPD] for other papers on 
the stochastic Navier-Stokes equation.) Therefore the regularity of the ini- 
tial velocity is important as well as the time-space regularity of the driving 
force. We shall provide existence and uniqueness results (global in time), 
as a generalization of the classical results in Hilbert spaces (recalled in the 
remark at the end of Section [5]) . 

The content of the paper is as follows. In the next section, we shall 
introduce the Navier-Stokes equation and the Besov spaces to work with. 
In Section O we shall deal with existence and uniqueness results of solutions 
to our problem on a small time interval, when the initial velocity belongs 
to a Besov space of negative order and the forcing term is integrable in 
time to some power with value into a Besov space of negative order in 
space (see Theorem 13. 1[) . In Section H] global (in time) existence results will 
be given; for this aim, first we shall split our problem into two auxiliary 
problems (considering the additive splitting of the velocity field u = x + y). 
The equation for the variable y will have small initial data and small forcing 
term; in Proposition 14.21 existence of a unique solution y will be proven. The 
problem for the variable x will be solved in Proposition 14.41 using a priori 
energy estimates, as in the classical case discussed e.g. in [Te79| . Theorem 
14.51 will concern the main result for the given problem in the variable u. 
Finally, we shall make some remarks on the case of a smooth forcing term in 
Section [5j The range of variability of the parameters involved in defining the 
functional spaces will be specified in Appendix A and B, providing examples 
of non classical solutions. 
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2 The Navier-Stokes equation 



We consider the Navier-Stokes evolution problem, i.e. the system of equa- 
tions governing the motion of an homogeneous incompressible viscous fluid 

8u 

— - vAu + (u-V)u+\7p = ip ^ 
V • u =0 

where the spatial domain is the torus T 2 = M 2 / (2-kZ) 2 . In other words, 
we consider our problem on the square [0, 2tt] 2 with periodic boundary con- 
ditions. We consider a finite time interval [0, T}. Here the unknowns are 
the velocity vector field u = (ui(t, £), u%{t, £)) and the pressure field of 
the fluid p = p(t,£) (for (i, £) G [0, T) x T 2 ); A is the Laplacian opera- 
tor A = 8 2 /<9£ 2 + <9 2 /<9£|, V is the gradient operator V = (8/8^1,8/8^2) 
and • is the scalar product in M 2 . (p is a given external force and v > is the 
(constant) viscosity coefficient. V • u = is the incompressibility condition. 

The basic functional spaces to set the problem are the generalized Sobolev 
spaces (s G M, I < p < 00) defined to be "regular" subspaces of the space 
T>' of periodic divergence-free vector distributions. To define them, we pro- 
ceed in the following way. Take any u G T>' . Then since V • u = on 
T 2 , there exists a periodic scalar distribution ij) on T 2 , called the stream 
function, such that 

u = V X V = (-8^/8^,8^/8^). (2) 

Decomposing tp in Fourier series with respect to the complete orthonormal 
system in L 2 (T 2 ) given by {^e lk <} k&2 

. ^ e ik < 

m = E 1^ 

fcez 2 

by ^) we get that u has the following Fourier series representation 

u(0 = E u k e k(Z), Uk G C, u k = -n_ fe , (3) 

where et(i) = ~2^\kl e%k ^- Note that {efc} fcgZ 2 is a complete orthonormal 
system of the eigenfunctions (with corresponding eigenvalues |/c| 2 ) of the 
operator -A in [L^ iv (T 2 )] 2 = {u G [L 2 (T 2 )] 2 : f j2 u(Q = 0, V • u = 0, 
with the normal component of u being periodic on 8T 2 }, k 1 - = (— A^, 
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\k\ = y/k% + k 2 and Zl = {k G Z 2 : |&| ^ 0}. 

Each e k is a periodic divergence-free C°°-vector function (i.e. e k 6D). The 
convergence of the series ([3]) depends on the regularity of the vector function 
u, and can be used to define Sobolev spaces as in the following definition. 
For any s G R, 1 < p < oo, we define 

H s p = { u G V : u = Y, u ke k and ^ u k \k\ s e k G [L p (T 2 )] 2 } (4) 

k&l fcezg 

i/p is a Banach space with norm = || Ylk u k\k\ s &k\\L p - 

Thus the unknown velocity will be considered as a function of the time 

variable taking values in some space 

u(t) = ^2u k (t)e k 
k 

Let P be the following projection from the space of periodic distributions 
onto the space of periodic divergence-free distributions: Pu = Y2k( u ' e -k)ek- 
Here (•,•) is the ([C°° (T 2 )] 2 )' - [C°°(T 2 )] 2 duality bracket. Applying the 
projection P to the first equation (P) we get rid of the pressure term (because 
(Vp, e_fc) = — (p, V-e_fc) = for any k), obtaining the following formulation 
of our problem 

u'(t)+Au(t) + B(u(t)) = f(t), te(0,T], (5) 

as an equality in the distributional sense, with some initial condition u(0) 
assigned. We have taken the viscosity v = 1, without loss of generality. 
A is the stokes operator A = —A, B is the quadratic operator B(u) = 
B(u,u) defined by the bilinear operator B(u,v) = P[{u ■ V)u]. Notice that 
B(u,v) is equal to P[V • (u <8> v)} because of the divergence-free condition. 
f = P<P- 

We are interested in the evolution problem ([5]) for initial data and forcing 
term not too regular (in space variables). For this purpose, we define a 
scale of spaces B pq consisting of the Besov spaces of periodic divergence- 
free vector fields. They can be introduced as real interpolation spaces (see, 
e.g. [BL] Theorem 6.4.5): 

B s ; q = {H s p \H s p *) e , q s* = (1-9)8! + 982, 0<6<1 (6) 

for s* G R, 1 < p, q < oo. 

In particular, B? 2 = H( and B% 2 = [L^ iv (T 2 )] 2 . 
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Since {efc} fcGZ 2 are the eigenvectors of the Stokes operator, with cor- 
responding eigenvalues = \k\ 2 , then we can naturally extend A to the 
whole space T>' by the following formula: Au = Uk\k\ 2 ek for any u £ T>' 
and by means of this representation it is straightforward to show that the 
Stokes operator A is a linear operator in Hp with domain Hp +2 ; moreover 
it is a bijective unitary operator from Hp~ onto Hp for any index s G R, 
1 < p < oo, and also from B^ 2 onto Bp q (for any s € R, 1 < p, q < oo) (see, 
e.g., [E] Theorem 1.1.6 for getting linear operators in interpolation spaces). 
In particular the inverse operator A" 1 is a linear bounded operator in each 
space Bp g . Finally, the operator A generates an analytic semigroup in each 
Bpq- (We refer, e.g., to [Te83] for classical analysis of the Navier-Stokes 
equation on the torus.) 

The properties of the Stokes operator are the basis to analyze equation (|5j) 
as a perturbation of the linear Stokes problem by the nonlinear operator B. 
To this end, the key point is to estimate the operator B in Besov spaces, as 
done, e.g., in (Ch96'J. 

Remark . We believe that all what follows holds if the spatial domain is any 
bidimensional smooth bounded domain D0, but we postpone this analysis 
to a future work. However, we point out that only in the space periodic case 
the expression of the Gibbs measure constructed by means of the enstrophy 
(given in [ARH-KJ) is an invariant measure for a stochastic Navier-Stokes 
evolution problem. This depends on the fact that this Gibbs measure is in- 
variant for the Euler flow both on the torus (see }ARH-K] ) or in any smooth 
bounded domain (see [AH-K] ), but the Euler and Navier-Stokes boundary 
conditions are the same only for D = T 2 . (For an overview on this subject 
see e.g. fAFa].) For this reason and because of the interest described in the 
introduction, we present our problem considering the spatial domain T 2 . We 
point out that our technique does not apply in the case the spatial domain is 
the whole space, because the inverse of the Stokes operator is not a bounded 
operator in Bp q (R. 2 ). Hence, the techniques are quite different from the case 
of the spatial domain D = R 2 . Moreover, also the assumptions for getting 
existence and uniqueness results are different from the case ID = R 2 ; let us 
consider [BGJ, which is as far as we are aware of the only paper including 
a forcing term in the analysis of very rough solutions of the Navier-Stokes 

1 For instance, there exists a complete orthonormal system {e n } nS N of eigenvectors of 
the Stokes operator in the space of square integrable divergence-free vector fields satisfying 
homogeneous Dirichlet boundary condition, with associated eigenvalues < Ai < X2 < • • ., 
An ~ n as n — > 00 (see, e.g., |Te79| ). This allows to analyze the linear Stokes operator as 
presented above for the case on the torus. 
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problem. Biagioni and Gramchev assume F G Lj oc (R + ; L q (M. 2 )),q G [g,2[, 
or F G L^R^L 1 ^ 2 ) n L°°(R 2 )) (setting F = V ± ■ f, for / the forcing 
term in ([5])). 

Other results for D = M 3 are presented in |CP| . □ 

3 Local existence and uniqueness 

In this section we prove a result on uniqueness and local existence for equa- 
tion ([5]) in Besov spaces. We use a technique from |B91j based on a theorem 
of local diffeomorphism of [VF| (§1.1). 

Let u be a vector valued function defined on the time interval [0, T] for 
T < oo or on [0, oo) if T = oo. Keeping in mind the notation in equation 
([!]) , consider the mapping $ 9 u i ► (u'+Au+B(u),u(0)) G £ -^F, for some 
Banach spaces The Frechet derivative d a <& at the point 

a = is given by the linear operator do<3? : £ — ► J-,u i— > (u' + An, u(0)) . Un- 
der the assumption that the quadratic operator B from £ to JF 1 is bounded^ 
and that the linear operator do& is an isomorphism, Vishik and Fursikov 
[VF] observe that <J> is analytic in a neighbourhood of G £ and locally $ 
has an inverse operator which is analytic in the neighbourhood j6(<7oo) 
of $(0) G T of radius (Too = (4||(d $)" 1 || 2 ||-B||) _ • This provides global 
existence of a unique solution u G £ for small forcing term and small initial 
data in T . This solution has analytical dependence on the initial value and 
forcing term. We will use this technique in Proposition 14.21 
Moreover, starting from this approach, [B91] shows a result of local existence 
for any data in J- by showing that the Vishik-Fursikov procedure works also 
for T < oo and by finding an estimate on ot- (The sub- index T reminds 
that we work on the time interval [0, T\. The important point is that also 
the norms in £ and in T will depend on T.) Roughly speaking, equation ([5]) 
is seen as a perturbation by a small nonlinear term B of the linear equation 
(which is well posed in the Hadamard's sense, see Proposition 13. 2p . The 
"smallness" of B is obtained choosing the time interval small enough. More 

2 In practice, these spaces will be constructed from L a (0, T; B") and Bp q spaces. From 
now on, Bp q denotes the space for vector valued functions defined on the torus, unless 
otherwise specified. By the way, we remark that in the notation for the time-integrability, 
the sup-index has been used, whereas in the notation for the space-integrability the mf- 
index has been used. 

3 The operator B : £ — > T 1 is bounded if 

||B|| := sup ||B(tt)||^i < oo. 
IMI £ <i 
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precisely, let the time variable vary in the finite interval [0, T]; assume that 
the quadratic term B : £ —* T 1 has a norm bounded by 

II-BII < Ci T £ 

for some positive constants C\ and e; then the solution u to the Navier- 
Stokes equation (JS} exists locally in time for any forcing term in T 1 and 
initial velocity in J-°. Indeed, given (/, uq) G J 7 , let < T < T be such that 

»(/■«•)»' T * < 4||(4,*)-|I'C. (?) 

Then (/, ito) € B(<7™), the .F-ball of radius cr^; according to Vishik and Fur- 
sikov's result quoted above, we conclude that there exists a unique solution 
defined on the time interval [0,T]. 

These are quite general results. The crucial point is the choice of the 
spaces £ and T . In |B91| regular Hilbert spaces are considered: £ = 
{u G L 2 (0,T;tff) : u' G L 2 (0,T;H$)} and T = L 2 (0,T;H$) x Hf. Here 
we consider £ = {u G Z/(0, T; B~^ +2 ) : u' G Z/(0, T; B~^)} and J" = 

— s+2— - 

L r (0,T;B~q ) x 5 pr r with 1 < p, q,r < oo and s£i The norms are 
defined as 

II (/, «o) IU = (Jo" ll/(<)llB-.^) 1/r + Kll R -» +2 -| 

^p? B pr r 

\He = (Jo \Wt)\\ r B s +2 dt + ; T n«'(t)||» B _,dt) 1/r 

q D p q 

We are especially interested in the case — s + 2 — - < 0. 

In this framework, the main result in this section is as follows. 

Theorem 3.1 Let 

£ = {ue 17(0, T; B--+2) : G Z7(0 5 T; B~«)} 

T 1 = L r (0,T; Bp* ) , T a = B~l +2 ~^ , F = F 1 x F° 
with 1 < p, q,r < oo, s£i. 
T/ten 

TTie map g?q^ : £ 3 ii i— ► (u' + Au,u(0)) £ J- is well defined and 
continuous. Moreover, the operator do<& is an isomorphism of Banach 
spaces £ and T . 
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ii) If the parameters r,p,q,s satisfy the following conditions 

r<q (8) 

*-l<! (9) 

1-8 <j (10) 

*+f+'<3 (11) 

2 < s + | + 7 (12) 
3<* + f + ^, (13) 

then the map £ 3 u i— > B{u) 6 J 71 is well defined and bounded. There- 
fore : £ — > J- is analytic. Moreover, there exist positive constants 
C\ and e such that 

\\B\\ < Ci T e . 

Hi) Under the same assumptions as in ii), for any forcing term f G 

L r (0, T;B~y) and initial data uq G B pr r there exists a time inter- 
val [0, T] C [0, T] and a unique solution u to equation dSJ) with initial 
velocity uq, defined on the time interval [0, T] with T > satisfying 
((7j), and 

u G L r (0,T;B-*+ 2 ) nC([0,T];B~: +2 - 2 -) 
u' €L r (0,T;B~°) 

The solution u depends analytically on the data uq and f. 
The properties i) of the linear operator do& follow from a general result: 



Proposition 3.2 Let T G (0, oo] ; 1 < p, q, r < oo and s G R. 

— s+2— 2 

For any / G L r (0, T; B~r) and uq G -B pr r ; i/iere exists a unique u £ 
W 1 » r (0,T) = {u£ L r (0,T;B p ° +2 ) : u' G L r (0, T; B~^)} such that 

(u\t)+Au(t) = f(t), i€(0,T] 
| u(0) = it 

Moreover, the functions u', u depend continuously on the data f and uq, that 
is there exists a positive constant c such that 

(/^(ii«(t)ii^-. +a +ii«'(oii^-.)*) 1/r < (cjo\\fm r B -»dt) 1/r Hu \\ z 
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Finally, the space W 1,r (0, T) is continuously embedded in the space C&([0, T];B P r 
that is there exists a positive constant c such that 



s+2- 



_ 2 



'C b ([0,T];B p r r ) 



< C ||«||wi,r( 0) D (14) 



and therefore the initial condition makes sense. 
All the constants^ depend only on r,p,q,s. 

For the proof, see, e.g., [B95| Proposition 4.1 (based on |D V] ) . The assump- 
tions on the linear operator A and on the space are fulfilled; namely, 
the properties of the Stokes operator recalled in Section 2 and the fact that 
B~q (1 < p, q < oo, s € K) is a UMD Banach space (i.e. has the Uncondi- 
tional Martingale Difference property). The last part of Proposition 13.21 is 
obtained by interpolation, bearing in mind the interpolation result for Besov 

spaces B pr r = (B~g, B'~|+ 2 ) 1 _i r for 1 < r < oo, 1 < p, q < oo, s € M. 

This corresponds to the first part i) stated in Theorem 13. 1[ Assuming i) 
and ii), then part Hi) is proven as described at the beginning of this section. 
Hence, the proof of Theorem 13.11 is complete as soon as we prove ii). This 
is given by the following result for the nonlinear equation. 

Proposition 3.3 Suppose the real numbers s and p,q,r E (1, oo) satisfy 
the conditions [M 13\) of Theorem \3.1\ Then there exist constants e > and 
C\ > such that for all T > 

\\B( U (t))\\ r B -sdt\ < dT s \\u\\ 2 wl , r{0>T) , ueW^&T). 

Proof. First, let us show that there exists a pair (a, b) of real numbers such 
that 



s < a,b (15) 
a,b<l (16) 



r 

2 
P 

2 



a + b=j + l-s (17) 

|[(o + 6) + 2(s-2 + 2)]<l (18) 

Before discussing these inequalities, let us remark that ([9|) and (fT7|) imply 
that 

a + b>0 (19) 



4 We make the convention to denote different constants by the same symbol c, unless 
we want to mark them for further reference. 
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We begin from (|12p . written as 



2 - 2 - s < 2 



v 



which grants that there are solutions to (|15p -(fT6 
Secondly, (fT8l) is equivalent with 



a + b<A-2s-i- (H8 



Since by (jlip (written as | + 1 — s < 4 — 2s — 2 ) any solution to (|17p satisfies 
(jl8fj) and hence (fTB|) . we only need to show that the system (fTI^ - ffTUjl - ffTTP 
has at least one solution. We look for a solution such that a = b. (Anyway, 
it is not difficult to see that a = b is not the only possible solution.) For this 
it is enough that 

Z r 6 ^ p + 2 2 ^ p 

The second of this inequalities reads 

which is (fTO)) . 
The first one reads 

which is equivalent to f)13[) . Thus, system (|15M18p has at least a solution. 
Define 

a = §[a + s- 2 + |] 
/3=5[fc + fl _2+a] 

By (I15p we have a, (3 > and by (|18p we have a + (3 < 1. In particular 
a,/? < 1. 

We are now ready to finish the proof of Proposition 13.31 First we estimate 
the bilinear operator by means of Bony's paraproducts techniques, as given 
in [Ch96], Corollary 1.3.1. Because a and b satisfy (fT6l) - (fT7l) - (fT9j) . we have 

\\B(u(t))\\ B -s = ||V-[«(t)®«(t)]||^. 

< ||«(t)®«(t)|| B -*fi (20) 

< c H<)l|B»,||w(t)||flg a 

Secondly, we use well known results on Besov spaces as interpolation spaces 
(see, e.g., [BL]) to get 



-Dp g 
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I II /II II 1-/3 II 11/3 
\u\\ R 6 < C U 2 W \\ „ 



-s+2 



with the interpolation parameters a, (5 defined above. Here and in the fol- 
lowing, c denotes different constants. 

We use these inequalities to continue the estimate of the quadratic operator 
from the last line of (|2(jp : 

\\B(u(t))\\ s < c\\u(t)\\ 2 ~ a J *\Ht)\\l + J s+2 

(21) 

< c\\u(t)f- a J 2 \\u(t)\\ a+ ^ +2 , for r<q 

Dp r 

Since a + (3 < 1, then 

(tf\\B(u(t))\\ r B -sdt) 1/r < c\\u\t a - P s+2 2 (J T \Ht)\\^f + \ r dt) 1/r 

< c\\u\\ 2 -^ s+2 2 T^-^(J^\\umi-^dt) {a+Pyr 

C([0,T];B pr + ') * Pq 

< C IImII 2 , , T (l-«-/3)/r 
— L II "II W l < r (Q,T) 

where we have used Holder's inequality for the time integral in the second 
line and the embedding (|14ft of Proposition 13.21 in the third line. ■ 

We point out that conditions ([9]) and (jllH are against each other; namely, 
rewritten down for the regularity value of the initial data, they are 

-S + 2-2>-2_(2_ 1 \ d _ s + 2 -2 > 2_ 1 

r r \p I r p 

It follows that 

-. + 2-2 >max{(-f -(§-!)), (|-l)} 

By the computations (|40p in Appendix A, we get that — s + 2 — - > — lat 
least. This imposes a restriction on the admissible initial velocity to solve 
equation © locally in time. 



4 Global existence and uniqueness 

We want to show that the local solution constructed in the previous section 
exists on the whole time interval [0,T]. To prove this, we split our problem 
into two subproblems, considering two auxiliary variables x and y such that 



11 



u = x + y. 

Following |GPj . let decompose the data as 



u = x + y and / = g + h 

The problem for the variable y will have small forcing term h and small 
initial data yo. Time-global existence and uniqueness will be proved by 
means of Vishik and Fursikov's technique. 

The problem for the variable x will have more regular data: initial data 
xq G V and force g £ L 2 (0,T; H^ 1 )- Time-global existence and uniqueness 
will be proved by means of an a priori estimate of the energy. 

By the very definition of Besov spaces ©, the space of periodic diver- 
gence-free smooth functions T> is dense in any Bp q . Therefore, we have the 
following Lemma for the splitting of the data. 

Lemma 4.1 Let f G L r (0,T; B° q ) with 1 < r,p, q < oo and a G R. Then 
for anye>0 there exist functions g £ G C°°([0, T);V) and h £ G L r (0, T; B° q ) 
such that ||/i £ ||l''(o,T;_b^ ij ) < £ an d f = 9 s + h £ ■ 

The case of constant (in time) functions gives the splitting for the initial 
data. 

We proceed now in this way. The two subproblems read 

' x'{t) + Ax(t) + B(x(t),x(t)) 

+B(x(t),y(t))+B(y(t),x(t))=g(t), * 6(0,2] (22) 

x(0) = Xq 

and 

/ y'{t) + Ay(t) + B(y(t), y(t)) = h(t), t G (0, T] 

\ 1/(0) = y (23) 

For the latter one, we can choose h and yo sufficiently small (written as ^SC 1 
below), in order to have the following result. 

Proposition 4.2 Let the assumptions [Ml3\) of Theorem \3.1\ be satisfied. 

— s+2— — 

Then, given any yo G B pr r , h G L r (0,T;B~ q i ) with \\yo\\ _ s+2 -2 1, 

Bp r 

ll^llL r (o t-b~ s ) < ^ 1> ^ ere exists a unique solution y to equation (j23[) on the 
time interval [0,2], such that 

y G U-(0,T;B^ 2 )nC([0,T};B; S r +2 -^), 
y'eL r (0,T;B-*). 
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Proof. Global existence for small initial data and small forcing term is 
obtained as described at the beginning of section [3] by means of Vishik and 
Fursikov's technique. Therefore, parts i) and ii) of Theorem 13. II entails this 
Proposition. ■ 

Concerning equation (|22p . we would like to show existence and unique- 
ness on the whole time interval [0, T]. We already know that there exist 

Te (0,T] and a function x(=u-y) e L r (0, T; B^ +2 )nC([0,T\; Bpr +2 ~ 7 ), 
x' € L r (0, T; B~q). To show that x is indeed a solution to ([22]) is enough to 
show that all the terms in (|22p make sense. This is done, for the nonlinear 
terms, analogously as in (|20l) - (|2T1) . For instance 

l \B(x,y)\\ B -s < c\\x\\ B aJy\\ Bb _ < INI 1 "^ \\x\\ a +2 ht'^ WvllL-** 



d r q d a ~ Ji ' r D p q 

JJ p q D p q 



2 2 2 2 
Condition (111]) says that — I h s < 3 (i.e. — I 1 < — s + 2 and 

p r p r 

- - 1< -s + 2 - -). Hence B^ +2 C B^q' 1 and B P r +2 ^ C Bj~ l and 
therefore 

3T € (0, T] and x such that 

x e L r {Q,T-Bjf r ~ l ) nC([0,T]; J B| r " 1 ), (25) 

i'a r (o,T;B;, +r3 ) 

Now we look for a priori estimates. By [GPj (Lemma 1.1), we have the 
following energy estimate 

2 2 

Lemma 4.3 Lei 2<p<oo,2<g<oo, — | 1>0. Then, for any 

p q 

e > there exists a constant c £ > smc/i that 



T 



(B(x(t)),y(t))dt 



T i-T 



< £ I \\x(t)f Hl dt+c t 

2 ; Q 



/ Mt)\\ 2 H 4y(t)\\ q 2+2 x cft. 



A precisation on the proof is required. In fact, Gallagher and Planchon in 
[GP] work in the whole space. Anyway, the technique used by them for 
the spatial domain D is M 2 , works also in our case D = T 2 . Indeed, let 
u = Yl,k u k e k with efc = efc(£) defined for £ G M 2 ; the Littlewood-Paley 
decomposition (see, e.g., |Ch98j and references therein) gives 

A m -u = ^2 u k e k 

2 m <|fc|<2 m+1 
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where A m means the convolution with a function ip m whose Fourier trans- 
form tp m has support in {£ € M 2 : 2 m < |£| < 2 m+1 }. Hence, the proof of 
[GP] based on Bony's paraproduct and Bernstein's inequality, is valid also 
if we deal with the Besov spaces 

B s pq (T 2 ) = {u = J2 w u k e k 6 V : £ meN (2 ma || A m u|| iji(T2) ) 9 < 00} 

and the L p (T 2 )-norms appear at the place of the L p (R 2 )-norms of |GPj . 
Moreover, the norms ||Vu|| i2 ( T 2) and ||it||iyi(x 2 ) are equivalent. Indeed, T 2 
is a bounded domain and (u, 1} = for u 6 V. 

We now look for a priori estimates for the unknown x. Let us multiply 
both sides of the first equation (|22l) by x and integrate in space and in time. 
Two terms vanish, namely (B(y(t), x(t)), x(t)) = and (B(x(t), x(t)), x(t)) = 
0, see e.g. [Te79] . Moreover (B(x(t),y(t)),x(t)) = -{B(x(t),x(t)),y(t)) and 
(Ax(t),x(t)) = ||x(i)|| 2 l . We then have 

±MT)\\ 2 h°+[ \\x(t)f Hi dt = ±\\x \\ 2 H o+ [ (B(x(t)),y(t))dt+ f (g(t),x(t)) 
2 Jo 2 2 Jo Jo 

< „ iNHLo+c / IkWIIIslll/WII'a+a,!* 

(by Lemma IQt 2 JO 2 B p q 

pq 



+ \ I Mt)\\ 2 H idt + c [ \\g(t)\\ 2 idt (26) 
io 2 Jo 2 

In this way, the required bounds are obtained by means of Gronwall's lemma, 
as soon as we can find p > 2, q > 2 with | + | — 1 > 0, such that 



T 



\y(t)\\ q a+a-i* < °°> 



_B?. « 

P 9 



where y is the solution to problem (|23p . Proposition 14.21 provides y € 
L r (0,T;Bp g s + 2 ). Thus we need to show that 



L r (0,T;B^ +2 )CL^0,T;Blp *) (27) 



for some p>2,q>2 with 4 + | — 1 > 0. 

5 See, e.g., [Ch98j for the definition and properties of Bony's paraproduct and [N] for 
Bernstein's inequality. 
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Further conditions on the parameters r, p, s are required in order that 
CUD holds. 

First, there is the embedding L r (0,T) C 1,9(0, T) if 

r > q (28) 

since the time interval is hnite. 



2,2 



On the other hand, the space embedding £~ s+2 C Bi s 9 ~ holds if 



pq — ^pq 

s + 2-|>|-l (29) 



1 < Q < 9 < oo, 1 < p < p < oo 

(see [BLj Theorem 6.5.1). 

We recall assumption (jlip of Theorem I3.lt 

s + - + - < 3 

p r 

The choice g = r, in order to satisfy (|28p . makes that the conditions (jlip 
and (|29p become identical (to be precise, there is the difference < or <, 
which makes (jlip a slightly stronger than (|29p ). Thus we have 

\\y\\ i+i-i < c \\y\\ B - s + 2 (30) 
up f 

Summing up, choosing q = r and p = p in (|26p and assuming ([8l ll3p with 



with the additional conditions r>2, - + - — 1 > 0, the proper estimates 



2,2 
p "r" r 

follow. More precisely, from (|26p we first have that 
lk(T) ||^o < ||ar ||#o + c / ||a;(t)||^o \\y(t)\\ r s+2 dt + c [ \\g(t)\\] U dt 

n 2 n 2 J n 2 Hpq J Q H 2 

Gronwall's lemma gives 

JJP IWHflS - C (H X 0||//0, l|y|lL-(0,T;B^+ 2 )' II^II^^.H-l)) < OO 

The last result in conjunction with ([26]) gives 



L 



10 

Finally 



T 

x(t)\\ 2 H i dt < CdlxollflO, l|y|lL--(0,T;B-=+ 2 )' M (0,T;H^)) < 00 



Il°°(0,T;H°) + \\ x \\l 2 {0,T;H^ - C ( 1 1 ^0 1 1 //0 , ||y |lir( ,r;Bp| +2 ) ' llfl'IIi, 2 (0,T;H^ 1 )) 

(31) 
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Use now the embedding theorem in Besov spaces 



a _ cr cr-l+Z cr-1+2 

H 2 = B 22 ^ B p2 P ^ B PP " 

(p>2) (p>2) 

so to get from (f3Tj) that 

l|X|l ^ (0 , T;B |- ) + ^ll L2{ 0,T;4) " C(l|XOll ^°' Hf IIl»(0,T;V)) 

2 2_j 

By (complex) interpolation (that is [L 2 ,L°°] 1 _2 = Z7 and [Bp° g ,Bp q ]j_2 = 

r ?' 

2 + 2_l 

Sp 9 r for 2 < r < oo), we obtain that 

ll ^ l iO o (0jT;B |- 1 ) + l|x|l ^(0 J T; B |^- 1 ) - CillX4H ^ llyll ^(°^ +2 )' I^W^)) 

(32) 

Comparing (|32p with (|25p . we get that the solution a; exists on the whole 
time interval [0, T\. We have therefore proven the following result. 

Proposition 4.4 Let the assumptions 151 113)) of Theorem \3.1\ be satisfied 
and moreover assume r > 2 and - + ~ — 1 > 0. Taen given any xq 6 P and 

a G i 2 (0, T; H^ 1 ), there exists a unique solution x to equation (j22|) on iae 
izme interval [0, T] sac/i i/ia£ 

x G 77(0, T; n C([0, r]; fl| r _1 ) 

3;'eL r (0,T;5;, +r3 ) 

Remark . Notice that more regularity on g does not improve the regularity of 
x, because of the presence of y (which has a role similar to an external force 
in equation (|22p ). For this reason, we assume g G L 2 (0, T; H^ 1 ) instead of 
the other possible choice g G C°°([0, T]; T>). On the other hand, the initial 
data xo is chosen very smooth in order to consider without problems the 
continuity in time in the next results. □ 

We combine Proposition 14.21 and Proposition 14.41 and, bearing in mind 
the embedding used to show (|25p (that is to show that x is less regular than 
y), we get that 



x + y e L r (Q,T; Bp°q~ *) D C([0, T];B£ r *) 

2,2 o 



u 



' = x' + y' £ L r (0, T; Bp° q 
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This implies that there exists a function u, given by u = x + y, with the 
regularity specified in ([33]) . This is indeed a solution to equation ©. In 
fact, bearing in mind equations (|22p and (|23p . we notice that the function 
u = x + y solves the equation © , thanks to the fact that the nonlinearity 
B(u(t)) is well defined; and this is so, because Chemin's result to estimate 
the quadratic term garantees that this exists if u belongs to some Besov 

space of positive index and from (|33|) we have that u(t) G Bp q r with 
| + ^ — 1 > (for a.e. t). Hence, u is the sought global solution. 

We sum up all the results proven so far and state our main theorem. 

Theorem 4.5 For any forcing term f G L r (0, T; B~t) and initial velocity 

— s+2— 2 

uq £ B pr r with 

s£l, 1 < q < oo 
2 < r < q 
Z + 2 - 1 > 

p r 

< s - 1< 2 
P P 

2<S +|+f< 3 

3<S+2+4 
p r 

i/jere exists a unique solution u to equation (UJ) on i/ie iime interval [0, T] 
u G L r (0,T;s|^ _1 ) nCaO.T];^" 1 ) 



Moreover, there exists a (strictly) positive T <T such that the above solution 

belongs to L 2 (0, T; B~^ +2 ) nC([0, T]; S pr r ). Hence strong continuity for 
t — > holds. 

Remark. In Appendix B, examples fulfilling all these assumptions will be 
given. The assumption r > 2 imposes that —s + 2 — 2 > —I as shown by 
(jllj) in Appendix A. □ 



Proof. What remains to be proven is the uniqueness result. Let us denote 
by U the set of functions u satisfying the conditions (|33p . Consider two 
solutions u, u £ U and denote by 5 the difference u — u. Then 5 €LU and it 
satisfies the equation 

U'{t) + A5{t) + B(u(t), 5(t)) + B(S(t),u(t)) = 0, t G (0, T] 

1<J(0)=0 
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This is a linear equation in 8. We analyze this equation as a linear Stokes 
problem with a (linear) perturbation term. If the perturbation B(u,-) + 
B(-,u) is good enough (mainly, small for small time T, so that this gives a 
small perturbation of the well-posed linear parabolic equation), then there 
exists a unique solution. This will hold on a small time interval; but since 
8 = is a solution, then we get that the unique solution is the zero one on 
a small time interval. Starting again from the zero value, we can proceed in 
the same way to cover the whole time interval. 

We want to analyze the perturbation B(u, 8) + B(S, u). We define the oper- 
ators 

F u S:=-B(u,S) T u 8 := -B(8,u) 

and the space 

S := { 5 € L r (0, T; : 5' € L r (0, T; b|, + '~ 3 ) } 

equipped with the norm || <5|| 5 = (/ T ||<S(t)|| r 2 + 2 .dt +fi \\8'{t)f 2 + 2 .Jt) 1/r . 

Dp" 1 " r 1 RP r 

D p q -Op q 

It is enough to consider the case with u, since the same works for u, because 
of the symmetry of Chemin's estimates in the two arguments. We are going 
to show that 

T u :S^L r (0,T;B}^) 

and 

IWH | +M < C U (T) \\8\\s (35) 

L-{Q,T;B p q r ) 

with C U {T) -» as T -> 0. 

This in nothing but an application of Chemin's estimates. In fact 



R p ' 



D P 9 


-2-2 

r 




2,2 
B p r 


,IH 


1 2-1 
R p 


B p 


,IH 


R P 



<c||u|| 2^2_,||(5|| 2_, for l + l-l >0 



< c [I ?x|| 2 , 2_, IIJII 2_, for r < o 

nP T r 
C P 9 

Therefore, integrating in time, we get 



\B(u,8)\\ 



2,2 



+ £-3. 



Lr(0,T;B p q r ) 



— C II U II 2 + 2-1 ||<5|| 2_! 

L-(0,T;B p q r ) C([0,T];S p p r ) 



<C 2 \\u\\ 2 + 2_, \\8\\ s . (36) 

Lr(0,T;B p q r ) 
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In the last step, we have used the interpolation result (as (|14p ) which allows 



to dominate the norm in C([0, T]; Bp r ) by the norm in S 

\\u\\ 2_! < C \\u\\s- (37) 

C b ([0,T];Bg r ) 



r-.g 



Finally, ([35} holds with 

C U {T) = C2IHI 2 + 2_i • 

Lr(0,T;B£ q r ) 

We remark that C U (T) -> as T -> 0. 
We shall show that the problem 

f <S'(t) + A$(t) = -r u 5(i) - lV(i), t E (0, T], 

\<5(0) = 

has a unique solution 5 E 5 on a small time interval, using a contraction 
theorem. Let us denote by T the mapping giving the solution to the Stokes 
problem 

(v'{t)+Av(t) =g 
\v(0) = 

We know from Proposition 13.21 that T is an isomorphism from 
L r (0,T;B^^~ 3 ) onto S and 

IMs < C 3 \\g\\ 2 +2 _ 3 (39) 

f(0,T;B/ a r ) 

for some positive constant C3. Therefore, by ([35]) and (|39|) we get that 

||T o (T u 6 + fuS)\\ s < C 3 (C U (T) + Cu(T)) \\S\\ S 

This shows that the mapping T o (T u + Tq) : S — > S 1 is a contraction as soon 
as we work on the time interval [0, T] C [0, T] with T > chosen in such a 
way that 

CsCO^ + C^T)) < 1 

Hence there exists a unique solution 5 £ S to equation (f34"j) on the time 
interval [0, T]. This must coincide with the zero function: 6(t) = for all 

te[0,T}. 

Since the constants providing the contraction mapping depend only on the 
norms of u and u (because problem (|34p is linear), we start again from 
S(T) = and we get the same result on the time interval [T, 2T] and so on 
to conclude the proof in a finite number of analogous steps. ■ 



19 



5 Remarks on the case where the forcing term / 
is smooth 

In the previous section, the constraint r > 2 has appeared. This implies 

-s+2-- 

that the initial data uq is assumed to belong to the Besov space B pr r 
with — s + 2 — | > — | (see the comment at the end of section [3] and (jHJ in 
Appendix A). This restriction comes from the use of the energy estimate of 
Lemma [4.31 We want now to show that for 1 < r < 2, when the forcing term 
/ is smooth enough ( say, / € L 2 (0, T; H 2 ) ), a classical energy estimate can 
be used instead of Lemma 14,31 Therefore, assuming conditions (l8l fT3j) and 
some more regularity on the forcing term, we prove global existence also 
when 1 < r < 2. Notice that, when the condition r > 2 is removed, the 
initial velocity regularity index — s + 2 — - can be very close to —1 (see 
also some examples of admissible values in Appendix B). This agrees with 
similar results obtained when there is no forcing term (see, e.g., |GP] for the 
problem in M 2 ). 

We go back to the statement of Proposition 14.21 assuming that h = (having 
taken g = /). Since y G L r (0, T; B~g +2 ), then there exists a t\ G (0, T] as 
close to as we want, such that 

When the condition 1 < r < 2 is added to ([8l ll3p . then 

-s + 2 > 

because at the end of section [3] we have shown that — s + 2 — ^ > — 1. 
Moreover, if the index regularity of the initial velocity is negative, i.e. — s + 
2 — £ < 0, then condition (fTT]) imposes p > 2. At the end of Appendix 

B, it will be shown that this implies B pr r % B® 2 - Therefore we are 
not dealing with the classical problem of initial velocity with finite energy. 
Summing up, if 1 < r < 2, -s + 2 - \ < and flj DSj) hold, then 

Since — s + 2 — 2 < < — s + 2, there exists 6 6 (0, 1) such that 

r° — (n~ s+2 ~'r n- s + 2 \ 

Therefore 

y(h) G B° p2 C B° 22 = Hi 
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because the spatial domain is bounded. (In fact, in the same way we can 
show that y is infinitely smooth in space and time on the time interval (0, T] , 
because there is no forcing term for y. But we do not need this result.) 
Therefore on the time interval [ii,T], the classical Hilbert-space theory can 
be applied to get existence and uniqueness results (see, e.g., |Te79j ): 

y€LL 2 (t 1} T;Hl)nC([ti,T\;H%) 

By interpolation between L 2 (ti,T; H),) and L°°(ti, T; H® ), we get 

y G L 4 (ti,T; H^ 2 ) 

We now assume / G L 2 (0, T;£?2 ). Since 

B° 2q CB~° if p > 2,-1 > -s- | 
L 2 (0,T) C L r (0,T) ifl<r<2 
Then in our setting we have that 

f GL 2 (0,T;B° 2q )cL r (0,T;B p °) 

Therefore local existence results are obtained by means of Proposition 13.31 
We only need an a priori estimate on the time interval [t\ , T] for the un- 
known x. This is easily obtained from the following classical estimate on 
the trilinear term 

(B(x),y) < \\x\\l4 ||Va;||i 2 ||y||i 4 by Holder inequality 

< c ||x[| 1/2 |la;|| rri |ly|| 1/2 by Sobolev embedding 

< c\\x\\ l jrl\\x\fjj 2 \\y\\ 1/2 by interpolation 

n 2 n 2 H 2 

< e \\x\\ 2 H i + c £ 11x11^0 112/11^1/2 by Young inequality 

Hence in the time interval [ti,T], the unknown x does not explode in the 
required norms and therefore there exists a unique x G L 2 (ti, T; H),) D 
C([t 1} T\;H$), x' G L 2 (t\, T; H^ 1 ). We remind that t x can be chosen close to 
as much as we want. Since existence on any small time interval [0, t{\ was 
already proven in section [31 this result implies the global existence. Finally, 
u G L 2 (h,T; Hi) n C([t%,T]; H%) for any < t% < T. And this solution u is 
unique. 

Remark. If uq G H® and / G L 2 (0,T; H^ 1 ), then a classical result grants 
that there exists a unique solution u G C([0,T\; H%) D L 2 (0,T; H^), u' G 
L 2 (0,T;H^) (see, e.g., [MS]). ' " □ 
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Appendix A Lower estimates on — s + 2 — ^ 

Let 1 < p, q, r < oo be given. Then 

' ] - 1 1 < p < 2, Vr 



max • 



(_2_2 + 1 2 
1. r P ' P J 

Hence, it easily follows that 



p>2, Kr<^ 



2 2 



- + 1 P > 2, r > 



inf max < - 

Kp<oo I 
KKoo 



2 1 + 1, I -l} 



P 



With some more (but elementary) work, we obtain 



mf max 

Kp<oo 
2<r<c» 



/_2 _2 +1 2_ 1 1 = _i 
[ r p ' p J 2 



(40) 



(41) 



Appendix B Admissible values for the parameters 
s,p,r and numerical examples 

We want to show that system (J9l ll3p has a non void set of solutions. These 
are the conditions appearing in Theorem 13. 1\ providing local existence. Two 
more conditions are required for the global existence of Theorem 14. 5( unless 
the forcing term is smooth enough (see Section [5]) . We start analyzing the 
less restrictive conditions (|9lfT3l). 

Set x = | and y = ~. Then the system of conditions is 

2 — s<x + y<3 — s 

3 — s < x + 2y 

s-Kx (42) 
1 — S < X 
< x, y < 2 

Because of the range of values specified in the last line, necessary conditions 
for the existence of a solution to (1421 ) are 

2 
3 
s 
1 



s < 4 and < 3 
s < 6 
1 < 2 
s < 2 
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Hence, the admissible values for the parameter s are — 1 < s < 3. Moreover 
x < 3 — s (from the first and the last line in (|42p ) and x > s — 1 (that is the 
third line in (|42p ) imposes the further restriction: s < 2. 
Summing up, the admissible values for the parameter s are 

-1< s < 2 



We distinguish two cases. 

• - 1 < s < 1 

Since s — 1 < 1 — s, the third line in (j4*2j) can be neglected. Representing the 
remaining conditions (|4"2|) on the (x, y)-plane, it is easy to see that there exist 
solutions. If we are interested in the solutions satisfying also the condition 
— 1 < — s + 2 — y < (for the regularity of the initial velocity), then s must 
be positive. We give examples of parameters satisfying the above conditions. 
Examples: 

for Pro. E31 s = ^ r=§ p = 12 - s + 2-2 = -| 
forTh.[43] s = -± r=f> p=§ - s + 2-2 = || 

We do not choose q since this is the less significant parameter to characterize 
a Besov space. Notice that 1 < r < 2 in the first case, providing global 
existence for "regular" forcing term. The second case concerns positive 
index regularity — s + 2 — - for the initial velocity. 

• 1 < s < 2 

Since s — 1 > 1 — s, the forth line in (j42]) can be neglected. Again the 
graphic representation shows that there are solutions. The condition — 1 < 
— s + 2 — y < (for the regularity of the initial velocity) requires x < 1 (i.e. 
p > 2). Notice that in this case both conditions y < 1 (i.e. r > 2) and 
— s + 2 — y < (i.e. — s + 2 — 2 < 0) can be fulfilled, providing parameters 
satisfying the assumptions of Theorem l4.5l with the initial velocity in a Besov 
space of negative order (not allowed in the previous case). 
Examples: 

for Pro. Ml s = %, r = f, p=f, - a + 2-2 = -g 
forTh.S3] S = if, r = f>, p = 4, - s + 2 -2 = -^ 

11 r _40 o _ s _L 2 -2-_J- 

s = |, r = 3, p=|, - s + 2-2 = o 
S = g, r = 21, p = 2, "5 + 2-2 = ^ 
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This is the only case providing global solutions with initial velocity uq in 
Besov space of negative order — s + 2 — £ < and force / € L r (0, T; B~^). 

Remark. We remark that -s + 2 - ^ < and s + | + £ < 3 imposes | < 1, 

that is p > 2. This implies that there the embedding B pr r C B 22 
never holds. Hence we really deal with a generalization of the classical 
result for uq £ -^2 2- An analogous statement holds for the forcing term: 
L r (0, T; £ L 2 (0, T; S^ 1 ), because £p 9 s £ S^ 1 for p > 2 and s > 1. □ 
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